A SENSITIVE VLA SEARCH FOR SMALL-SCALE GLYCINE EMISSION 

TOWARD OMC-1 


J. M. Hollis 1 , J. A. Pedelty 2 , L. E. Snyder 3 , P. R. Jewell 4 , F. J. Lovas 3 - 5 , Patrick Palmer 6 , 

and S.-Y. Liu 7 - 8 


1 Earth and Space Data Computing Division, Code 930, NASA Goddard Space Flight 
Center, Greenbelt, MD 20771 

2 Biospheric Sciences Branch, Code 923, NASA Goddard Space Flight 
Center, Greenbelt, MD 20771 

department of Astronomy, University of Illinois, 1002 W. Green Street, 

Urbana, IL 61801 

4 National Radio Astronomy Observatory, P. O. Box 2, Green Bank, WV 24944-0002 
5 Current address: National Institute of Standards & Technology, 100 Bureau Drive, 
Stop 8441, Gaithersburg, MD 20899-8441 
department of Astronomy and Astrophysics, University of Chicago, 5640 South Ellis 
Avenue, Chicago, IL 60637 

department of Astronomy, California Institute of Technology, Pasadena, CA 91 125 
8 Current address: Institute of Astronomy and Astrophysics, Academia Sinica, 

P.O. Box 23-141, Taipei 106, Taiwan, R.O.C. 



- 2 - 


ABSTRACT 

We have conducted a deep Q-band (X -7 mm) search with the Very Large Array (VLA) 
toward OMC -1 for the lowest energy conformation (conformer I) of glycine 
(NH2CH2COOH) in four rotational transitions: the 615 - 5 14 , 6 24 - 5 2 3, 7 17 - 6^, and 
7 q7 - 6q 6 • Our VLA observations sample the smallest-scale structures to date in the search 
for glycine toward OMC- 1 . No glycine emission features were detected. Thus if glycine 
exists in OMC- 1 , either it is below our detection limit, or it is more spatially extended than 
other large molecules in this source, or it is primarily in its high energy form (conformer II). 
Our VLA glycine fractional abundance limits in OMC -1 are comparable to those 
determined from previous IRAM 30 m measurements — somewhat better or worse 
depending on the specific source model — and the entire ~T primary beam of the VLA was 
searched while sensitive to an areal spatial scale -150 times smaller than the 24 " beam of 
the IRAM single-element telescope. In the course of this work, we detected and imaged the 
4 ] 4 - 3j 3 A and E transitions of methyl formate (HCOOCH3) and also the 2 0 2 -1 01 

transition of formic acid (HCOOH). Since formic acid is a possible precursor to glycine, 
our glycine limits and formic acid results provide a constraint on this potential formation 
chemistry route for glycine in OMC- 1 . 

Subject headings : ISM: abundances — ISM: clouds — 


ISM: individual (Orion Molecular Cloud) — 
ISM: molecules — radio lines: ISM 
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1. INTRODUCTION 

Recent laboratory experiments demonstrate the synthesis of amino acids, a basic constituent 
of life, when ice mixtures representative of interstellar grain mantles were subjected to UV 
radiation at temperatures below 15 K. Bernstein et al. (2002) produced amino acids in a water- 
rich ice environment, while similar results were achieved by Munzo Caro et al. (2002) in a 
water-poor ice environment. Sorrell (2001) suggests glycine could be produced from simple 
interstellar species that undergo radical-molecule interactions when excess reaction heat is 
available to overcome the activation barrier. Sorrell (2001) proposes the following reaction 
sequence may occur in grain mantle sites: 

CO + OH -4 COOH (1) 

CH 3 + COOH -4 CH3COOH (2) 

NH 2 + CH3COOH -4 NH 2 CH 2 COOH + H (3) 

Such interactions would require high concentrations of requisite radicals and molecules that 
may be present in such sites. Moreover, acetic acid (CH3COOH) has now been detected in two 
interstellar clouds: SgrB2(N-LMH) and W51 (Mehringer et al. 1997; Remijan et al. 2002). 
Another possible route for producing glycine is a combination of gas-grain chemistry involving 
formic acid (HCOOH), a known interstellar molecule, and the amino-alcohol (NH2CH2OH) 
which has been identified as a component of Comet Halley dust (Kissel & Krueger 1987). 
Protonated glycine is produced in an endothermic reaction (Chamley, Ehrenfreund, & Kuan 
2001) which would find favorable conditions in hot molecular cores: 

NH 2 CH 2 OH£ + HCOOH -> NH 2 CH 2 COOH 2 + H 2 0 (4) 

Additional support for the probable existence of interstellar glycine is provided by the recent 
detections of interstellar glycolaldehyde (CH2OHCHO) and ethylene glycol (HOCH2CH2OH) 
by Hollis et al.(2000, 2002) which emphasize that large biomolecules are common constituents 
of interstellar clouds. 
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The first glycine spectrum observed in the laboratory was a high energy conformer (later 
designated conformer II), and its assignment was reported contemporaneously by two 
independent groups (Brown et al. 1978; and Suenram & Lovas 1978). Later, Suenram & 
Lovas (1980) analyzed conformer I glycine. Lovas et al. (1995) showed that conformer I 
glycine is 705 cm' 1 (1012 K) lower in energy than conformer II glycine, and that the dipole 
moments are p a = 0.91 1 D and pb = 0.607 D for conformer I and p a = 5.372 D and pb = 
0.93 D for conformer n. (While at least six glycine conformers are possible in theory, only 
conformers I and II have ever been observed in the laboratory.) Since molecular line 
intensity is proportional to the square of the dipole moment, glycine a-type transitions are 
the more likely to be detected in interstellar clouds. Conformer n spectra dominate 
conformer I spectra in a laboratory absorption cell because of its larger dipole moment and 
the high sample cell wall temperature (>450 K). The expected line intensities of the two 
conformers equalize at ~285 K. Thus, if the temperature of an interstellar molecular source 
is < 100 K (e.g., OMC-1), then conformer I glycine a-type transitions would be the 
candidate lines most likely to be detected. 

Following the successful laboratory measurements, searches for interstellar glycine have 
been conducted by various groups (Brown et al. 1979; Hollis et al. 1980; Snyder et al. 

1983; Beralis et al. 1985; Guelin & Cemicharo 1989; Combes, Rieu, & Wlodarczak 1996; 
Snyder 1997; Ceccarelli et al. 2000). Combes et al. (1996) conducted the most sensitive 
previous search for glycine using the IRAM 30-m telescope at 3 mm with a 24" beam size; 
at the same time, they also searched at 1 and 2 mm with beam sizes of 13” and 19", 
respectively. We proposed to conduct a glycine search with the NRAO 8 Very Large Array 

9 The National Radio Astronomy Observatory is a facility of the National Science 
Foundation operated under cooperative agreement by Associated Universities, Inc. 


(VLA) for a number of reasons: First and foremost, interstellar glycine may have a spatial 
scale of one or two arc seconds; such a scale size can only be detected with a sensitive 
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interferometer because beam dilution would plague a single-element antenna. Second, 
telescope sensitivity is proportional to collecting area and the VLA has -20 times the 
collecting area of the IRAM 30-m telescope. Third, all VLA antennas have been recently 
equipped with receivers at Q-band (-7 mm) which is the only VLA operating band with 
enough tunable range (40-50 GHz) in which several likely transitions of glycine could be 
observed. Fourth, the VLA acts as a spatial filter which resolves out emission from small 
molecules that form in extended envelopes, thereby reducing the line confusion level to 
make it possible to detect larger molecules that are prone to survive in embedded compact 
sources. Fifth, different complex molecules often have very different spatial distributions; 
such spatial association can be used to help sort through the spectral line confusion. 

This work reports the results of a deep search of the low temperature environment of 
OMC-1 with the VLA for four conformer I glycine a-type transitions in the Q-band range. 
This is the first Q-band glycine search to be conducted with an interferometer. 

2. OBSERVATIONS 

Simultaneous Q-band observations of glycine and silicon monoxide or methyl formate 
and formic acid were conducted toward OMC-1 with the VLA in 2001 September (DnC 
configuration) and 2001 December (D configuration) at epochs and scheduled observing 
time specified in Table 1 cols. (1) and (2), respectively. The J2000 phase center for all 
observations was a = 5 h 35 m 14 s .25, S = -5°22'35".5, which is midway between the Orion 
hot core and the OMC-1 formic acid positions given in Liu et al. (2002). Moreover, the 
phase center position is within ~1" of the pointing position used by Combes et al. (1996). 
The VLA primary beam is « 45/v(GHz) arcminutes - thus, the nominal primary beam for 
our observations was -T. The VLA correlator was operated in the Four IF normal mode, 
which yields 32 spectral channels with a spacing of 195.313 kHz (-1.36 km s* 1 ) for both 
left and right circular polarizations. An assumed LSR source velocity of +8.0 km s' 1 was 
centered in the resulting 6.25 MHz (-43.5 km s' 1 ) bandpass. Observed molecules and 
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bandpass center frequencies are shown in Table 1 cols. (3) and (4). Table 2 lists the 
glycine transition parameters: the rotational quantum numbers (col. 1); the transition rest 
frequency (col. 2); the upper energy level E u (col. 3) in Kelvins; and the transition line 
strength S (col. 4). 

Quasars 0607-085 and 3C84 were used as the phase and bandpass calibrators, 
respectively. The absolute flux density scale was determined from observations of 3C286 
(Sep 16 through Dec 20) or 3C48 (Dec 23 and 26) whose 43 GHz flux density is 1.44 Jy 
or 0.52 Jy, respectively, as given by the SETJY routine in AEPS. The latest antenna gain 
curves were applied using the FIXQGAIN procedure (S. Meyers, private communication). 
The bootstrapped flux densities of 0607-085 and 3C84 at each observing epoch (Table 1 
cols. 5 and 6, respectively) show satisfactory repeatability and are in good agreement with 
NRAO monitoring results. However, on the dates where 0607-085 and 3C84 fluxes 
appeared systematically low or high (see Table 1), a median flux value of 1.63 Jy is adopted 
for 0607-085. 

The antenna gains were derived from 0607-085 observations at 10-30 min intervals and 
applied to the data using standard AIPS routines. Phase self-calibration on vibrationally 
excited SiO maser emission (Table 1) was used to derive fast timescale (5 s) atmospheric 
phase variations and the solutions were applied to the corresponding glycine line data. This 
correction was not applicable for the methyl formate and formic acid data. The extensive set 
of SiO observations will be the subject of a future paper. 

Continuum subtraction was performed using the UVLIN routine. Continuum subtracted 
data cubes of total intensity emission were created by averaging the right and left circularly 
polarized visibilities. The rms noise levels in each plane of a continuum subtracted cube 
range from 15 to 50% above the theoretical noise limit, which is reasonable for a low 
declination source (maximum elevation was ~ 50°). The naturally weighted synthesized 
beams resulting from DnC and D configurations were ~1.9"x 1.6"and ~2.9"xl.9", 
respectively. 
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3. RESULTS 

During the 2001 September observations, we devoted each track to one transition of 
glycine (as shown in Table 1) except the 7 0 7 - 6 0 6 transition which was observed twice 

because it has the largest line strength (see Table 2, col. 4). Based on BIMA Array formic 
acid observations (Liu et al. 2002), our search probed regions of peak formic acid emission 
which are found predominantly along the OMC-1 compact ridge at +8.0 km s* 1 LSR 
velocity. In 2001 December, we observed the 7 07 - 6 0 6 glycine transition in three tracks 

with a four channel frequency shift. Table 2 col. (5) summarizes the emission limits 
provided by our glycine search. Two remaining tracks in 2001 December were used to 
detect the 2 02 - lot transition of formic acid and the 4 J4 - 3 j 3 transition of methyl formate. 

Formic acid was chosen since both it and glycine share the structural carboxyl group 
(COOH), which is seemingly rare in interstellar clouds. Because of this structural 
similarity, it may well be that both species share common formation pathways. Further, one 
would expect to find both species in spatial proximity, but not necessarily spatially 
coincident, if the formic acid is a precursor which is used up in the process of glycine 
formation. Methyl formate was chosen because it is a ubiquitous interstellar isomer of 
both acetic acid and glycolaldehyde, and its 8-atom size is more comparable to 10-atom 
glycine. Table 2 cols. (5), (6), and (7) summarize the formic acid and methyl formate 
emission features detected. 

Table 2 lists the beam-averaged column densities Nt (cols. 9 and 10) for formic acid and 
methyl formate and upper limits for glycine. These column densities were calculated from 
the following expression (Liu, Mehringer, & Snyder 2001) which assumes optically thin 
LTE conditions and does not account for possible beam dilution: 

N _ 2 17xl0 20 AI AV QROT exp(E u /TRo T ) 

e a e b v 3 sp 2 


( 5 ) 
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In equation (5), Nj is in cm' 2 ; v is the rest frequency in GHz; E u is the upper rotational 
level in K; S is the line strength; p is the dipole moment in debye (i.e., 0.91 1 D); AI is the 
peak line intensity in Jy beam -1 ; AV is the FWHM linewidth in km s -1 ; the product 0 a 0 b is 

the Gaussian FWHM dimensions of the synthesized beam in arcsec 2 (Table 2, col. 8); 
Qrot, the dimensionless rotational partition function, is 15.6xTrot> L7xTrqt, or 
12.45x Trot for glycine, formic acid, or methyl formate, respectively (e.g., Gordy and 
Cook 1984); and Trot, the rotational temperature, is assumed to be either 43 K, a value 
based on a rotational diagram method for formic acid (Liu et al. 2002), or 100 K, a value 
approximately representative of other molecular species toward the Orion compact ridge 
(e.g., Wright et al. 1996). 

Figures 1 and 2 show contour maps and spectra, respectively, of the A and E symmetry 
states of the 4j4 - 3 j 3 transition of methyl formate. Positions of continuum source I and 

the infrared BN object (e.g., Menten & Reid 1995) are provided in Figure 1 and subsequent 
spatial imaging figures for reference. Figures 1 and 2 of methyl formate demonstrate that 
the A and E symmetry lines possess nearly identical spectra and nearly the same spatial 
morphology as expected for high signal-to-noise data (cf. Figure 1 panels d and e of Minh 
et al. 1993). 

Figures 3 and 4 show a contour map and spectrum, respectively, of the 2 02 - Iqi 
transition of formic acid that results from a convolution with a 6" circular gaussian to 
enhance the extended structure of this weak emission. As a result, formic acid is 
convincingly detected at the ~20 mJy beam -1 level for this extended beam size. For a 
naturally weighted beam size, as shown in Table 2, formic acid has a peak intensity of ~5 
mJy beam -1 . Liu et al. (2002) have shown that methyl formate and formic acid in OMC-1 
are not spatially coincident. In Figure 5, methyl formate emission is between source I and 
the formic acid source. The position of continuum source I marks the centroid of SiO 
maser emission and is, therefore, the source of molecular outflow. Liu et al. (2002) suggest 
the formic acid source is located in a layer which delineates the shock interaction as the 
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outflow plows into the ambient quiescent gas. Thus, in this scenario, methyl formate would 
represent older, post-shock gas. 


4. DISCUSSION 

The detection and spatial imaging of the 4 14 - 3 j 3 A and E transitions of methyl formate 
(HCOOCH 3 ) and the 2 02 -loi transition of formic acid (HCOOH) in OMC-1 
demonstrate that these weak molecular emission lines originate in compact sources. Since 
formic acid is a possible precursor to glycine (as shown in equation 4), our glycine limits 
and formic acid results provide a constraint on this potential formation chemistry route for 
glycine in compact sources of emission toward OMC-1. For example, comparing column 
density results (Table 2 col. 9) for the same beam size (Table 2 col. 8 ), the abundance of 
glycine is less than or equal to the formic acid abundance in OMC-1. 

The most severe upper limit recently quoted (e.g. Ehrenfreund & Menten 2002) on the 
glycine-to-molecular-hydrogen fractional abundance ratio (Xoiy) toward interstellar sources 
is of order 10' 10 , which is based on observations with the IRAM 30m telescope (Combes et 
al. 1996). We will show that this value is optimistic. Any such quoted limit must be 
carefully qualified. First, conversion of a measured upper limit for a line flux to an upper 
limit for a beam-averaged column density depends on several assumptions. Typically, for 
upper limits, one assumes that the optical depth is small, LTE conditions at a known 
temperature prevail, and possible beam dilution is ignored. This is what we have done in 
equation (5) and Combes et al. (1996) have done in their work. Second, forming abundance 
ratios implicitly assumes that the two molecules are spatially coincident. Even with this 
assumption, the abundance ratios are model dependent unless both measurements were 
made with the same beam. For the same total number of molecules in the source, the beam- 
averaged column densities derived from observations that reach the same flux-sensitivity 
with two different beams will be different. If the source is small relative to the beam, there 
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will be beam-dilution, but, if the source is large relative to the beam, flux will be missed 
because the instrument is only sensitive to beam-sized spatial scales. 

Beam dilution must be addressed when comparing observations where telescope beams 
differ significantly (e.g., the VLA and IRAM 30m telescope) by computing a beam filling 
factor (B), a fraction <1 which is applied by division into the observed emission intensity 
(equivalently the column density as per equation 5 which assumes B=l). For a circular 
gaussian telescope beam size (0^) centered on the peak of a circular gaussian source of size 
( 0 S ), the resulting simple convolution yields a beam filling factor (B) given by (e.g., see 



For an infinitely extended source, B approaches unity. If the source size is very small 
compared to the telescope beam, B approaches zero. A source size must be assumed, 
estimated or modeled a priori so that the actual value of B can be determined for each 
telescope beam involved. Thus, in the following, we compare VLA and IRAM 30m Xoiy 
ratio limits that result from taking into account source size. 

Our VLA experiment is the first sensitive search for glycine in OMC-1 sources that have 
angular sizes on the order of -2". We consider three possible glycine source sizes in 
OMC-1 when computing Xoiy ratios. Since glycine contains an amine group (NH2), we 
choose nitrogen bearing molecules that have been observed with an interferometer to obtain 
physically meaningful source sizes. From Wright et al. (1996, figure 1), we select ethyl 
cyanide (CH3CH2CN), a nine-atom molecule that rivals the size of glycine; from Wilson et 
al. (2000, figure 2), we select ammonia (NH3), a well-known hot core constituent of OMC- 
1; and, finally, we assume a true source size equal to VLA beam size. Table 3 summarizes 
the apparent size (col. 1), the imaging beam size (col. 2), and the true size (col. 3) for these 
three selections. True sizes (s) were deconvolved from the observed apparent sizes (obs) in 
the circular beam (b) approximation (e.g., see equation 8.15 of Rohlfs & Wilson 2000): 
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e obs = e s + 0 b (7) 

To estimate the molecular hydrogen column densities needed to form Xdy ratios, we 
follow Wright et al. (1996) who use observations of optically thin 3 mm dust emission in 
OMC-1. Molecular hydrogen column density determinations from dust emission are not 
straightforward and have an order of magnitude uncertainty (Wright et al. 1996). For the 
OMC-1 compact ridge which was observed with a beam size of 7" x 5", Wright et al. 

(1996, table 4) derived an observed molecular hydrogen beam-averaged column density of ~ 
2 x 10 23 cm' 2 . 

VLA 3a upper limits to the glycine beam-averaged column density toward the compact 
ridge in OMC-1 (Table 2 cols. 9 and 10) assume a beam filling factor of unity and are a 
function of the synthesized beam size used, the rotational temperatures assumed, and the 3o 
noise levels estimated in the final images. The 7 07 - 6 q 6 transition was observed with more 

integration time than other glycine transitions, and therefore provides the best constraint on 
glycine total column density for a naturally weighted VLA beam size of 2.3" x 1.7" (see 
Table 2). The rotational temperature of glycine is problematical in OMC-1 where 
characteristic temperatures typically range from 50 to 100 K (e.g., see Combes et al. 1996). 
However, our VLA observations would be more sensitive to hotter compact core emission, 
and therefore we use the hotter Table 2 column density 3a limit (i.e., < 4.64 x 10 15 cm' 2 ) as 
characteristic of the transition of glycine for the 2.3" x 1.7" VLA beam case. The IRAM 
30m telescope was used at 3 mm (FWHM = 24") to obtain ala limit for a beam-averaged 
glycine column density of < 5 x 10 13 cm' 2 assuming a temperature of 100 K (Combes et al. 
1996; N.-Q. Rieu, private communication). For comparison of results between the two 
telescopes, we use the 3a glycine beam-averaged column density upper limits of <4.64 x 
10 15 cm' 2 and < 1.5 x 10 14 cm -2 for the VLA and IRAM 30m telescope, respectively. 

When forming Xciy from column density measurements or upper limits obtained with 
different telescope beams (Q), one must multiply the beam-averaged column density ratio 
by the inverse of the ratio of beam filling factors: 
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x _ N T( G] y) b(Q(h 2 » 

Gly N t (H 2 ) B(Q(Gly)) 

As per equation (8), beam filling factor ratios and Xoiy ratio limits are shown in Table 3 
cols. (4) and (5), respectively, for the three source models considered for each telescope 
(i.e., either the VLA or the IRAM 30m as shown in col. 6). 

Table 3 shows that for decreasing source sizes, the Xoiy ratio limits range from < 1.6 x 
10' 8 to < 4.7 x 10‘ 9 for the VLA while the same limits range from < 4.6 x 10' 9 to < 1.1 x 10- 
8 for the IRAM 30m. Alternatively stated: For a source size like that of ethyl cyanide, the 
IRAM 30m telescope produces an Xdy ratio limit ~3.5 times more severe than the VLA; on 
the other hand, for source sizes like that of ammonia and the VLA beam itself, the VLA 
produces Xoiy ratio limits of ~ 1.5 and ~2.3 times, respectively, more severe than the IRAM 
30m. 

It is emphasized that Table 3 limits are based on conformer I glycine observations only 
and at temperatures < 100 K. Under more extreme hot core conditions (e.g., 300 to 600 K), 
the prevalent form of glycine may well be conformer II (see discussion in § 1). 

5. SUMMARY 

We conducted a deep Q-band search with the VLA toward OMC-1 for conformer I 
glycine in four rotational transitions: the 655 - 5^, 6 2 4 - 5 2 3, 7j7 - 615, and 7 q 7 - 606- 
No glycine emission features were detected, suggesting that either conformer I glycine is 
below our detection limit, or is more spatially extended than other large molecules in OMC- 
1. Alternatively, glycine may be primarily in its high energy form (conformer II). For 
future OMC-1 searches, we found that the four glycine Q-band transitions we attempted to 
observe appear clear of contaminating line emission from other species. We carefully 
discussed the necessity to account for the sizes of the beams used in all of the observations 
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— both glycine and molecular hydrogen. Our measurements were used to determine upper 
limits on the total column density of glycine toward OMC-1. These results yield Xciy ratio 
limits ranging from of < 1.6 x 10' 8 to < 4.7 x 10' 9 for source sizes of 12.1" x 5.9" to 2.3" x 
1.7", respectively. We conclude that the VLA conformer I Xciy ratio limits in OMC-1 are 
comparable to the previous IRAM result -- somewhat better or worse depending on the 
assumed size of the glycine source — and in none of the cases have Xoiy ratio limits for 
either telescope attained a value as low as lO 10 . Furthermore, the VLA results apply to the 
entire ~T primary beam while being sensitive to an areal spatial scale -150 times smaller 
than the 24" beam of the IRAM single-element telescope. We were able to detect and 
spatially image the 4 j 4 - 3^ A and E transitions of methyl formate and also the 2 02 - lot 

transition of formic acid. These positive results demonstrate the viability of the VLA as a 
sensitive search instrument for weak molecular lines of biological importance. 

We appreciate the VLA scheduling efforts of Barry Clark on our behalf, the suggestion 
of an anonymous reviewer that we compare VLA and IRAM results for a number of source 
sizes, and clarifying information provided by Nguyen-Quang Rieu regarding previous 
glycine limit results with the IRAM 30m telescope. J.M.H. and J.A.P. received support 
from NASA RTOP 344-02-03-01. L.E.S. and F.J.L. received support from the Laboratory 
of Astronomical Imaging at the University of Illinois, and NSF grant AST 99-81363. P.P. 
received support from NASA grant NAG5-8708. S.-Y. Liu received support from NSF 
grant AST 99-81546. 



TABLE 1 


Summary of Observational Parameters 



Observing 

1st Line 

2nd Line 

0607-085 

3C84 

Date 

Time 

Frequency 

Frequency 

Flux 

Flux 

(2001) 

(1) 

(hrs) 

(2) 

(MHz) 

(3) 

(MHz) 

(4) 

(j y )a,b 

(5) 

(Jy) a ,b 

(6) 



NH 2 CH 2 C00HC 

SiO (J = l-0)d 



Sep 16 

5.5 

43183.651 

43122.080 (v=l) 

1.55(2) 

6.87(16) 

18 

7.0 

43753.841 

43423.858 (v=0) 

1.88(11) 

8.52(50) 

27 

5.0 

42985.223 

43122.080 (v=l) 

1.64(3) 

7.31(12) 

28 

6.0 

42637.789 

42820.587 (v=2) 

1.62(3) 

7.30(15) 

29 

6.5 

43753.841 

43423.858 (v=0) 

1.64(2) 

7.36(12) 

Dec 07 

5.0 

43754.621 

43423.858 (v=0) 

1.48(2) 

6.99(18) 

10 

10.0 

43754.621 

43423.858 (v=0) 

1.66(2) 

9.93(17) 

20 

5.0 

43754.621 

43423.858 (v=0) 

1.06(2) 

4.92(15) 



HCOOCH 3 e 

HCOOHe 



23 

3.5 

45396.636 

44911.75 

... 

... 

26 

4.0 

45396.636 

44911.75 

0.95(2) 

4.05(14) 


a Flux density is the average at the two observing frequencies. A 1.63 Jy flux density for 0607-085 was 
used in calibration procedures on Sep 18, Dec 20, 23, and 26 (see § 2). 

^Uncertainties in parentheses refer to the least significant digit and are la values. 

September glycine observational frequencies are line rest frequencies as shown in Table 2. December 
glycine observational frequencies are shifted higher approximately four channel spacings (see § 3). 
dRest frequencies from Lovas (1992). 

e Rest frequencies shown in Table 2. Methyl formate is approximately centered between the rest 
frequencies for the A and E symmetry states. 






TABLE 2 


Summary of the Spectral Line Search Toward OMC-1 


NtxIO-1 4 N t x10 ' 14 


Transition 

J fc-K+ -J K-K+ 

(1) 

Frequency 

(MHz) 

(2) a 

E u 

(K) 

(3) 

s 

(4) 

Al(peak) AV 

(mJybeanr 1 ) (kms* 1 ) 
(5)^ b (6)*> b 

VlsR 
(km s* 1 ) 
(7) a 

0 a x6 b 
C X ") 
(8) 

Trot =43 K 

(cmr 2 ) 

(9)3, c 

Trot=100K 

(cm -2 ) 

(10) a > c 

Glycine**: 










615- 5 ]4 

42637.789(6) 

7.574 

5.774 

<0.9 

3.9 


1.9x1. 6 

<35.4 

<114.0 

6 24 - 5 2 3 

42985.223(30) 

8.420 

5.350 

<0.9 

3.9 


1.9x1. 6 

<38.0 

<121.0 

7l7~6l6 

43183.651(25) 

8.700 

6.818 

<0.9 

3.9 


1.9x1. 7 

<27.9 

< 88.1 

7 07 ~ 6 06 

43753.841(28) 

8.657 

6.851 

<0.6 

3.9 


2.3x1. 7 

<14.7 

< 46.4 





<3.0 

3.9 


6.0x6.0 

<8.0 

< 25.2 

Formic Acid e : 










2 02 -1 01 

44911.75(5) 

3.233 

2.000 

5-0(3) 

3.9(7) 

8.0(7) 

3.0x1. 6 

13.1(32) 






20.0(10) 

3.9(7) 

8.0(7) 

6.0x6.0 

7.0(16) 


Methyl Formate*: 










4 14~ 3 13 e 

45395.775(16) 

6.146 

9.947 

85.8(4) 

3.9(7) 

7.7(7) 

2.9x1. 5 

274(51) 

894(166) 

4 14 -3i3 a 

45397.394(16) 

6.128 

9.947 

84.4(4) 

3.9(7) 

7.7(7) 

2.9x1. 5 

269(50) 

879(163) 


Uncertainties in parentheses refer to the least significant digit(s) and are 2 a values for col. (2) and la estimates 
for cols, (5), (6), (7), (9) and (10). 

b G!ycine limits in cols. (5), (9), and (10) are 3a estimates and glycine linewidths assumed in col (6) are based on the 
nominal values observed for formic acid and methyl formate. 

c Abundance values or limits were computed using equation (6) with Trot = 43 K (col. 9) or 100 K (col. 10) as justified in § 3. 
^Dipole moment p a = 0.91 1 D (Lovas et al. 1995); rest frequencies calculated in this work. 
e Dipole moment p a = 1 .391 D (Kim, Keller, & Gwinn 1962); rest frequencies (Willmot et al. 1980). 
f Dipole moment p a = 1.63 D (Curl 1959); rest frequencies (Oesterling et al. 1999). 





TABLE 3 

Glycine-to-Molecular-Hydrogen Abundance Ratio 3a Limits for Three Source Models 


Apparent 

Imaging 

True 

B(Q(H 2 )) 

B(£2(Gly)) 



Source 

Beam 

Source 

XGly 

Telescope 

(" x ") 

(" x ") 

(" x ") 




(1) 

(2) 

(3) 

(4) a 

(5) 

(6) 

14 x 7.7 

7x5 

12.1 x 5.9 

6.1445 

<4.6 x 10-9 

IRAM 

(CH 3 CH 2 CN map) 


0.6995 

<1.6 x 10-8 

VIA 

4.9 x 1.9 

1.3 x 0.9 

4.7 x 1.7 

13.2146 

<9.9 x IO- 9 

IRAM 

(NH3 map) 



0.2825 

<6.6 x 10-9 

VLA 

N/A 

N/A 

2.3 x 1.7 

14.9077 

<1.1x10-8 

IRAM 



(assumed) 

0.2010 

<4.7 x IQ 9 

VLA 


a The ratio of the beam filling factor of Wright et al. (1996) -- who used a 7" x 5" beam to obtain an 
observed Nt(H 2) = 2 x 10 23 cm' 2 - to the beam filling factor of either the 24" circular beam of the IRAM 
30m or the 2.3" x 1 .7" gaussian VLA beam used to observe glycine (see § 4 text and equation 8) 
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FIGURE CAPTIONS 

Fig. 1 — Contour maps of the A and E symmetry states of the 4 j 4 - 3j3 transition of 
methyl formate toward OMC-1. The maps were made directly from continuum subtracted 
u,v visibilities; data from spectral channels 10 through 13(18 through 21) were Hanning 
smoothed and summed to produce the E symmetry state (A symmetry state) contour plots 
in panel a (panel b). Comparison of panels a and b verify expectations that transitions of a 
given molecule with the same excitation parameters (i.e., line strength, energy level, dipole 
moment, etc.) should be morphologically identical. The 2.90" x 1.52" naturally weighted 
beam size is shown in the lower left of each panel. Contour intervals are -3.5, 3.5, 7, 10.5, 
14, 17.5, 21, 24.5, 28, 35, 42, 49, and 56 mJy beam -1 with peak fluxes of 60.2 (58.7) mJy 
beam-1 in panel a (panel b). In both panels, the 3.5 mJy beam' 1 contour level corresponds 
to a 3 ct detection level. 

Fig. 2 — Spectrum of the 4 14 - 3 j 3 transition of methyl formate showing A and E 
symmetry states toward OMC-1 at 195 kHz channel spacing. Spectral channels 4 through 
28 are shown. Transition quantum numbers are shown in the upper right. The frequency 
axis and the LSR velocity axis reflect the LSR velocity of -8.0 km s' 1 assumed for this 
source position at a rest frequency of 45396.636 MHz (i.e., approximately midway between 
the A and E features). This spectrum was taken at the most intense position shown in 
Figure 1 with abeam size of 2.90" x 1.52" . 
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Fig. 3 — Contour map of the 2 02 - Iqi transition of formic acid toward OMC-1. A map 
was made directly from continuum subtracted u,v visibilities and convolved with a 6" x 6" 
circular gaussian, whose size is shown in the lower left. The convolution was necessary to 
recover flux from this obviously extended source. Data from spectral channels 15 through 
17 were Hanning smoothed and summed to produce the contour plot. Contour intervals are 
-6, -3, 3, 6, 9, and 12 mJy beam' 1 with a peak flux of 13.6 mJy beam' 1 . The 3 mJy beam -1 
contour level corresponds to a la detection level. 

Fig. 4- Spectrum of the 2 02 - Iqi transition of formic acid toward OMC-1 at 195 kHz 
channel spacing. Spectral channels 4 through 28 are shown. Transition quantum numbers 
are shown in the upper right. The frequency axis and the LSR velocity axis reflect the LSR 
velocity of ~8.0 km s _1 assumed for this source position at a line rest frequency of 
4491 1.75 MHz. This spectrum was taken at the most intense position shown in Figure 3 
with a 6.0” circular synthesized beam. 

Fig. 5 - Contour map of the average of A and E symmetry states for the 4 )4 - 3 13 

transition of methyl formate toward OMC-1 superimposed on a gray scale intensity map of 
the 2 q 2 - loi transition of formic acid. See Figure 1 caption which describes the methyl 

formate map production and similarly Figure 3 caption regarding formic acid map 
production. Contour intervals for the methyl formate are -4, 4, 8, 12, 16, 20, 24, 28, 32, 40, 
48, 56, and 63 mJy beam' 1 with a peak intensity of 63.8 mJy beam' 1 . The 2.90" x 1.52" 
naturally weighted beam size for the methyl formate is shown in the lower left. 
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